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ABSTRACT 

Context. We have re-analyzed the X-ray spectra of the gravitational lensed high-redshift BAL QSO APM 08279+5255, observed 

with the XMM-Newton and Chandra observatories. Previous studies (Hasinger et al. 2002; Chartas et al. 2002) detected unusual, 

highly-ionized iron absorption features, but differed in their interpretation of these features, regarding the kinematical and ionization 

structure. 

Aims. We seek one physical model that can be successfully applied to both observations. 

Methods. For the first time we have performed detailed photoionization modeling on the X-ray spectrum of APM 08279-1-5255. 

Results. The absorbing gas in APM 08279-1-5255 can be represented by a two-absorbers model with column densities A'h(I) « 7x10^^ 

cm"^, Njj(2) a 6 X 10^^ cm"^, and ionization parameters log^(l) » 1.5 and log ^(2) » 3, with one of them (the high-ionization 

component) outflowing at v » 0. 18(±0.01)c, carrying large amount of gas out of the system. We find that the Chandra spectrum of 

APM 08279+5255 requires the same Fe/0 ratio overabundance (previously) indicated by the XMM-Newton observation, showing that 

both absorber components underwent similar chemical evolution and/or have similar origin. 

Key words, galaxies: active - X-rays: galaxies - quasars: absorption lines - quasars: individual (APM 08279+5255) 



1. Introduction 

Broad absorption line (BAL) quasi-stellar objects (QSOs) are 
objects displaying in their spectra broad (FWHM a; 10000 km 
s"') absorption lines in the rest-frame ultraviolet (UV), origi- 
nated in outflows of matter from the c entral engine of QSOs 
JFoltz et al.lll990l: IWevmann et all 1 19911). Th e outflow velocity 
may reach up to 0.2c (e.g.. iFoltz et alj[l983l) . Determining the 
relationship between the material absorbing the X-rays and the 
one absorbing the UV radiation is key to our understanding of 
the geometry and the physical state of the medi um surround- 
ing the vicinity of supermassive black holes (e.g..'Math ur et alJ 
11995; Murray et al. 1995; Hamann 1998; Proga et al. 200^^ 

Before ChandrajXMM-Newton missions, detections of BAL 
quasa rs in X-ray were rare. Usually, these object are X-ray weak 
(e.g., [Gallagher et al.l [2006), sometimes interpreted as strong 
excess absorption. Chandra and XMM-Newton observations 
of BAL QSOs, have provided new constrai nts in the amount 
of absorption toward selected objects (e.g., Sabra & HamannI 
l200UIOshima et al.ll200U [Gallagher et al]l2002ll2006l) . indicat- 
ing large column densities Nh <: 10^^ cm~^. 

The BAL QS O APM 08279+5255 at redshift z = 3.91 
dlrwinet alJll998l) is one of the most luminous objects in the 
universe, further magnified by gravitational lensing by a fac- 
tor of ~ 50 - 100 (e.g., Ledoux et a l. 1998). It was detected 
with the Submilimiter Common-User Bolometric Array, imply- 
ing an apparent far -infrared luminosity greater than 5xlO'^L0 
JLewis et al.l Il998l) . The optical spectrum, obtained with the 
High Resolution E chelle Spectrometer at the Keck telescope 
(lEUison et al.lll999 ). along with a detailed study of the physical 
condit ions in the BAL flow of the QSO bv lSrianand & PetitieanI 
( I2OOOI) . allowed them to conclude that the corresponding gas 



stream, outflowing with velocities of up to 12000 km s ', is 
heavily structured and highly ionized. 

The quasa r APM 08279+5255 was observed twice with 
XMM-Newton dHasinger et al J2002l hereafter H02). In both ob- 
servations the quasar is observed clearly out to 12 keV, which 
corresponds to almost 60 keV in the rest frame. The most ap- 
parent feature in the XMM-Newton spectrum is an absorption- 
like feature around 1.55 keV (which they interpret as an ab- 
sorption edge), corresponding to ~ 7.7 keV in the rest frame 
of APM 08279+5255. The high-inferred iron abundance at 
the high redshift, corresponding to a young age of the uni- 
verse, is of great interest in the context of chemical enrich- 
ment models, and provides constraints on the early star for- 
mation history of the universe and on its cosmological pa- 
rameters (e.g.. Hamann & Ferlandl [19931 ; iHasinger et al] 120021 ; 
Komossa & Hasinger 2003) . 

The q uasar APM 08279 +5255 was also observed with 
Chandra dChartas et alj 1200 2'. hereafter C02). The Chandra 
spectrum shows a similar absorption feature as the XMM- 
Newton observation, but the feature led to a different interpreta- 
tion. In particular, C02 modeled the spectrum with two absorp- 
tion lines at 8.1 keV and 9.8 keV in the rest frame of the quasar, 
interpreted as Fe xxv K lines. If the Chandra interpretation of the 
dataisright, the te/A: velocity of the X-ray BALs is ~ 0.2c-0.4c. 
The presence of similar outflow velocitie s has been claimed 
in a few other AGN X-ray spectra (e.g ., [Pounds et alj 120031 ; 
IChartas et al]|2003l; [Pounds & Pagell2006l) . but alternative inter- 
pretations of the same spectra have been propose d, which do not 
requi re these relativistic outflow velocities (e.g.. [Kaspi & Behan 

[2OO6I) . 

J^ew models (e.g., [Elvisl [20001; [Proga et al.[ [20001; [Progal 
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[2007) predict that a large fraction of the accreted matter 
into the region of the compact object, is expelled out again 
in the form of high-velocity outflows. Broad UV absorption 
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features, with velocities ~ 0.05c - 0.1c, are associated with 
this material through acceleration mechanisms like accelera- 
tion b y gas pressure (e.g., We ymann et al. 1982 ; Begelman et al. 
1991 1 ) due to dust (e.g., IVoit et alj 119931; lYun & ScoviUg 
1995h and acceleration du e to radia tive pressure by spectral 
lines (e.g. , JDrew & Bokse nber3 [198 4; Shlosman ^t all 1 1 9851; 
Arav et al. 1994; de Kool & B egeknanI 119951; [Murray et alj 
1995; Proga et al. 2000) ob served in about 10% of luminous 
high-z quasars dLaor & Bran dt 2002), implying that these out- 
flows are an important component of the general picture of 
AGNs. Furthermore, these models also predict that in order for 
the UV material to reach such high velocities a shield made by 
a high column density of gas (Nh ~ 10^^ - 10^"* cm"^) must 
absorb radiation in the X-ray band, in this manner preventing 
the destruction of the UV material. The new generation of X- 
ray observatories, XMM-Newton and Chandra give us an unique 
opportunity to study in detail the X-ray component of the BAL 
QSOs, contributing in our knowledge about the dynamical and 
physical evolution of the center of these systems, and in the case 
of high-redshited BAL QSOs, the gas enrichment history in the 
early universe (Hamann et al. 2004). 

However, none of the previous studies include a self- 
consistent photoionization modeling of the X-ray spectrum of 
APM 08279-H5255. Any constraint on the ionization level of the 
absorbing gas and its connection with the kinematical proper- 
ties of the BAL outflow is of major interest to elucidate some 
of the greatest discrepancies between these two proposed mod- 
els. Furthermore, a model that can separate the differences be- 
tween observations (XMM-Newton vs Chandra), or unify both in 
a single frame, is highly desirable. We report a spectral analysis 
(made of these two observations separated by ~ 2 weeks in the 
rest frame) of the high-redshifted BAL QSO APM 08279+5255, 
and present one model that might reconcile both observations in 
the same physical context, shedding new light on the ionization 
degree, kinematics and evolution of this system, in connection 
with its cosmological consequences. We use the following cos- 
mological parameters: Ho 
Q^ = 0.7. 



70 km S-' Mpc-', Q.m = 0.3 and 



2. X-ray observations and Spectral analysis 

The quasar APM 08279+5255 was observed twice wifli XMM- 
Newton. The first X-ray observation was made on 2001 October 
30 with ~ 15 ks of exposure time (referred as XMMl in Table 2 
of H02). A significantly longer (~ 100 ks) observation was made 
in 2002 April 28-29 (referred as XMM2 in Table 2 of H02). 
Taking advantage of the improvement in calibration and effective 
areas used in XMM-Newton observations, we re-processed the 
primary event file and extracted the spectrum, using the most 
updated XMM-Newton Science Analysis System (SAS, version 
7.0.0), and the standard data pipeline processing. In our anal- 
ysis, we used the ~ 88.8 ks Chandra ACIS-I observation of 
APM 08279+5255 (see Table 1 of C02), and extracted the spec- 
trum, built ancillary and redistribution matrix files with CIACQ 
version 3.4. Details on calibration fluxes and count rates can be 
found in the respective works. 

We will present the spectral analysis in three steps: 1) By 
showing results from C02 (below in this section), 2) comparing 
them with H02 (below in this section), and 3) we introduce our 
own approach highlighting the differences between these sets of 
data {XMM-Newton vs Chandra, in Sect.[3]l. Then we proceed to 



' Chandra Interactive Analysis 
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of 



Observations. 



three more steps: 1) fitting the edge and line models; 2) fitting 
photoionization models to the XMM-Newton and Chandra data 
separately; and 3) fitting jointly in the context of one consistent 
model (in Sect.|4]i. 

In C02, the authors are able to fit an absorbed power-law 
with an intrinsic neutral absorber in the frame of the quasar, to 
the Chandra spectrum of APM 08279+5255. They find a best- 
fit photon index F x; 1.7 and the column density of the ab- 
sorber Nh ~ 6 X 10^^ cm"^. Significant residuals are found at 
~ 8 keV (rest-frame). This is the source of an important con- 
troversy. By the time C02 was under review, the authors had 
published H02, in which a good description of the spectrum is 
proposed by fitting a model that included an absorption edge 
(zphabs X PI X 1-Edge or Edge model) at ~ 8 keV to the 
XMM-Newton observation of this object. Despite the fact that 
Chandra data provides values for the parameters for the Edge 
model, its;^'^ is statistically poor (~ 1.36 for 107 degrees of free- 
dom [d.o.f]). Moreover, C02 tried to model this feature with a 
Single-Gaussian line model (the one at ~ 8.1 keV) and the fit 
gets a bit worse (xl ~ 1.41 for 106 d.o.f). They found a better 
fit with a Two-Gaussian lines model at ~ 8.1 keV and ~ 9.8 
keV (rest frame), which indeed significatively improved the de- 
scription of the data over the Edge model. A visual comparison 
between the two sets of data suggests that the absorption feature 
is present in both cases, but the line-shape of the feature appears 
to be weakened in the XMM-Newton observation. 

As noted by H02, a major improvement (in the description 
of the XMM-Newton spectrum of APM 08279+5255) is found 
by adding an absorption edge to the absorbed power-law (Fit 2 
in Table [TJ. Our best-fit model for the XMM-Newton data pro- 
vides an edge energy at x 7.7 keV (in the rest frame) and an 
optical depth t x 0.3. This was interpreted by the previous au- 
thors as a Fe edge, compatible with ionization potentials of iron 
from Fe xv to Fe xvm, implying a significant ionization of iron. 
Indeed, in that work the authors compared the absorption seen in 
the spectrum at low-energies and arrived at the conclusion that 
an overabundance of iron (~ 2 - 5) with respect to lower Z el- 
ements (O, Ne, Mg, Si and S) is necessary to account for this 
feature (A^pe ~ 1.5 x lO''' cm^^). 

In all our fits, we have included a neutral abso rption with 
the G alactic value fixed to A^h = 3.9 x 10^° cm"^ (IStark et al.1 
Il992l) . We begin our spectral analysis by fitting an absorbed 
power-law with an intrinsic neutral absorber in the frame of the 
quasar ( z^ 3.91). O ur XMM (EPIC pn+MOS) and XMM EPIC 
pn data ( Striider et al.ll200lb . are consistent with each other and 
with the Chandra fit (included also the H02 fits for comparison 
purposes). The photon index is F ~ 1 .7 - 2 and the column den- 
sity of the absorber is Nh ~ 6 x 10^^ cm"^. Table [T] quotes the 
parameter values of the model along with the goodness-of-fit in 
terms of a global x~ for the d.o.f used for each dataset (bins in- 
cluded in the range 0.2 - 1 ke V minus the number of free param- 
eters of the model), and the gZofofl/;^^^ -probability {Py^^i), which 
gives the probability of exceeding x^ for v degrees of freedom. 
Additionally, we introduce a criterion to assess the fits close to 
the absorption feature, thus we take the global parameters of the 
model and compute ;i'^ in the band (1-2) keV (observed frame), 
along with its corresponding Zocfl/;i'^ -probability. 

We include in this analysis, fits of the XMM-Newton data 
of the Single and Two-Gaussian models, and the results are 
presented in Table [T] for two sets of XMM-Newton data. One 
in which we fit simultaneously the EPIC pn, MOS 1 and MOS2 
(pn+MOS) data and another only taking the EPIC pn data 
(pn). In the case of the Single-Gaussian line model, both 
sets are compatible with each other and with the Chandra fit. 
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Table 1. Fit Results and comparison between sets of data. XMM-Newton vs Chandra observations of APM 08279+5255. 



Parameter 




H" 


C^ 


(pn+MOS)-^^ 


pn'' 


Ctf 


Fit 1: PI and Neutral 


Absorption at Source 








r 

Norm' 

Nh (xlO^-cm-^) 

;r'/(d.o.f) 




2 ()4+().uj 

6-92!°| 
481.8/(365) 


1 7O+0.0J 

^•' -0.06 

6.0!«| 
182.1/(109) 


I 99+0.02 

1 90+0:02 
^•^■-'-0.02 

^ 9^+0.21 

564.5/(478) 


2 05 +"■"■* 

1.25ti 

6.28!«i^ 

330/(290) 


J 79+0.03 
^•'J^-002 

n 1 Q+0.41 

229.7/(184) 


/"globalLFit 1] 




2 X 10-^ 


1.5 X 10-^ 


3 X 10-3 


0.05 


0.01 


;r^(l-2keV)/(d.o 


■f) 






134.7/(59) 


72.4/(35) 


85.7/(52) 


/'(l-2keV)[Fitl] 








2 X 10-' 


6 X 10-5 


2 X 10-3 


Fit 2: PI and Neutral 


Absorption at Source, and One Edge 






r 

Norm' 

Nh (xlO^-cm-2) 

£edge (keV) 

T 

x'/(d.o.f) 




2oi+u.(u 

734^1 
7.68t!« 

394.6/(362) 


7.68!;-^' 
0.37i« 
146/(107) 


J 95+0.02 

1 28"*^"""^ 

6^20ti 
77Q+0.16 

30+°-'^^ 
"•-""-o.oe 

506/(476) 


2.02+'"" 
1.31+""^ 

"•■'^-0.26 

7,53+0-21 
292/(288) 


1 77+0.02 
^■"-0.02 
1 09+"°' 

7.50+11 

7.70^1' 

190.1/(182) 


/"global [Fit 2] 




0.12 


6 X 10-3 


0.17 


0.44 


0.35 


X^{\-2k£Y)l(±o 


■f) 






74.5/(59) 


35.4/(35) 


53.7/(52) 


Al-2keV)[Fit2] 








0.09 


0.51 


0.46 


Fit 3: PI and Neutral 


Absorption at Source, and One Gaussian 






r 

Norm' 

Nh (xlO^-cm-^) 






1 7^+0.06 
^■'■-'-0.06 

6 4"^"-* 
8.05fi? 
<0.140 

r, -,,+0.06 


1 QQ+0.02 
^•^='-0.02 
J 25+0.01 

gQ4+o:23 

1 9 +0:04 
°-^^-0.03 
0.34+013 

i4+!^:ii 

"■1^-0.11 


2,04+ooj 
J 26+o:o2 

6 28+"'^^ 

7 97+0:03 

0.i5+"07 

n 1^+0:22 
0.16_ojo 


i.80!i!i;^ 

1.06+""' 

7 39+0:25 

g J J +0:01 

QQ9;8:8i 

17+8:85 


-Elincl 
0"lincl 

EW,i„e, (keV) 






X'/Cd.o.f) 






146.9/(106) 


531.3/(476) 


310.8/(288) 


198.4/(181) 


/"global [Fit 3] 






5.3 X 10-5 


0.04 


0.18 


0.35 


;r2(l-2keV)/(d.o 


■f) 






98.8/(59) 


53.2/(35) 


55.29/(52) 


/'(l-2keV)[Fit3] 








5 X 10-* 


0.03 


0.4 


Fit 4: PI and Neutral 


Absorption at Source, and Two Gaussian 






r 

Norm' 

A'/y (xl02-cm-2) 
£ii„ei (keV) 
o"iinci (keV) 






1 77+0.06 

8.05!^ 
<0.140 


1.98+""- 
1.28+^1 

g 05+0:02 

9 53 
80+ 
iq+o:io 

509.7/(474) 


2,03+003 
J 29+0:02 

6-49+^-i 
793+8:63 

0.09+""^ 

0.30+^'''^ 

9 52+ 

80+ 
A 21+0:13 

"■^1-0.11 

293.2/(286) 


J 79+0.02 

J 08+""^ 
74710:24 

g o9+o:oi 

0.06+°°' 

0.31+^-^^ 
9 72+8:05 

07o+8:m 

A 99+0:14 

178.9/(179) 


EWunei (keV) 
Eiinc2 (keV) 






Q 24+0.06 

g 7(j+0:20 

Q4j;8:l? 
A 40+0:17 

106.7/(103) 


a-iine2 (keV) 
EWii„,2 (keV) 

;r'/(d.o.fi 






/"global [Fit 4] 






0.41 


0.13 


0.39 


0.52 


;r'(l-2keV)/(d.o 


■f) 






75.3/(59) 


35.5/(35) 


44.8/(52) 


/'(l-2keV)[Fit4] 








0.08 


0.51 


0.79 



All ab sorp tion-line and edge pa rameters are computed for the rest frame. The error parameters are 90% confidence limits. (a) lHasinger et alJ 
1 I2OO2I) (b) IChartas et alj ( l2002h (c) The fits are carried out fitting simultaneously the pn+MOS data. 293 PI bins are from the pn data, and 188 
bins are from the MOS data (in the range 0.2-10 keV). (d) Only pn data, (e) Our reprocessed Chandra data, with 187 bins, (f) Power-law 
normalization, xlO-'' photons keV-'cm-"s-' at 1 keV in the observed frame. Pl= Power-law. 



and the three sets are in agreement with an absorption feature 
with equivalent width (EW) ~ 0.2 keV. For the Two-Gaussian 
lines model, the best-fit parameters lead to absorption lines at 
~ 8 keV and ~ 9.5 keV (rest frame) with EW ~ 0.30 keV 
and ~ 0.20 keV, respectively (all of them compatible with the 
Chandra measurements, within the errors). But, in this case, 
the;^'^-probabilities are (for the Chandra fits) higher than those 
for the Single-Gaussian line model, and the Edge model 
(with 4 parameters more), in both senses, global and locally. 
In principle, the fits of both models (Fit 2 and Fit 4) are for- 
mally acceptable for both sets of data. A more careful analysis, 
using physically-motivated arguments is required to elucidate 
the ambiguity of this issue. A temporal variability study may 



help to solve some aspect of the problem. We present a brief 
overview of the variability observed in the X-ray spectrum of 
APM 08279+5255, as seen from the point of view of the XMM- 
Newton and Chandra data in the next section. A more quantita- 
tive physical discussion is given in Sect. |6] 



3. Variability of the absorption structure at 
keV - 8 keV 



7.7 



For this analysis, we worked on two sets of XMM-Newton data; a 
16 ks exposure time observation of APM 08279+5255 (XMMl), 
and a 100 ks exposure time observation (XMM2), see Table 1 of 
H02. As a first step, we were interested in searching for variabil- 
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Fig. 1. Chandra X-ray spectrum of APM 08279+5255 with 
two models over-plotted: the Edge model (solid red line) and 
the 2-Absorption-lines model (solid blue line). In the right 
panel, we show a closer view of the residuals in the range 
1-2 keV. Solid line is the Edge model. Dashed line is the 
2-Absorption-lines model. 



ity within each observation. We looked at the count rate light 
curves for XMMl and XMM2. We see slight changes in the 
light curve, although almost always within ~ 10%, obtaining 
an average of ~ 0.17 cts s"' for XMMl and ~ 0.13 cts s"' for 
XMM2. For the rest of the analysis we will only use the data of 
the XMM2 observation of APM 08279+5255, since it collects 
one order of magnitude more photons with respect to XMMl, 
increasing notably the signal-to-noise ratio, and we refer to it as 
the XMM-Newton data. We also looked for evidence of vari- 
ability in the Chandra observation of APM 08279+5255, but no 
significant deviation from the average ~ 6.6 x 10"^ cts s"' was 
found, fully compatible with the count rate reported by C02. 

Here, we discuss in some detail the differences we are see- 
ing between the XMM-Newton and Chandra observations, which 
are separated by ~ 2 months in observed frame (~ 2 weeks 
in the frame of the quasar), in the context of best-fit models. 
To begin with, we will assume that indeed both sets of data 
are better represented by two different models (Edge model 
vs 2-Absorption-lines model), and study how the Chandra 
data (better represented by the 2-Absorption-lines model, 
see C02 for details) behaves with the best-fit model parame- 
ters extracted from fitting procedure, using XMM-Newton data 
(specifically, the Edge model parameters taken from Table 2 of 
H02), and vicerversa (i.e., how the XMM-Newton data behaves 
with the best-fit model parameters found in fits using Chandra 
data). 



First, we applied the simplest strategy: We use ;if -statistics 
to describe the goodness-of-fit of the data to the model. In Fig. 
[1] (left panel), we plot the Edge model, taking the intrinsic ab- 
sorption Nh, optical depth t and energy of the edge from Table 
2 of H02, leaving only the normalization and the photon index 
r free to vary (all these numbers are fully consistent with our 
own best-fit parameters, see Table [1]). The resulting j^'^ is 192 
for 185 d.o.f. We also plot the 2-Absorption-lines model, 
but this time we allow all parameters free to vary (all the num- 
bers fully consistent with the parameters reported by C02). This 
time x^ - 179 for 180 d.o.f. At this stage, we are unable to 
judge which model is better at describing the data, since it is 



^ And from XMM2, we only make use of the EPIC pn data, 
since it has a better calibration below 0.8 keV than the EPIC MOS 
(see http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0018-2-4.pdf I. 
No significant differences are seen between MOS and pn data, and the 
conclusions of the analysis can be easily applied to the MOS data as 
well. 



true that the 2-Absorption-lines model reports a smaller ;if^, 
but it also introduces more free parameters in the fitting proce- 
dure, which might produce the advantage. A simple hypotheses 
F-test is not permitted under this context, partly because these 
two are not nested models (a detailed discussi on about condi- 
tions under which the F-test is valid is given in iProtassov et al.l 
,2002) . Therefore, because the differences between models are 
actually local, we focused on the spectral discrepancies close to 
the strongest absorption feature at ^ 1.6 keV (observed frame), 
specifically the range 1-2 keV. We have adopted both mod- 
els using the global best-fit parameters, but evaluate x^ locally 
from 1 to 2 keV. The results are: ;tf^[Edge](i_2keV) = 55 and 
;t'^[2Lines](i_2 keV) = 45 for 52 d.o.f. First, it is easy to see that 
the differences between globals;^'^ (mostly) come from this spec- 
tral region. The largest discrepancies between models are con- 
centrated in the region around to 1 .6 ke V. We show these diff'er- 
ences in Fig. [T] (right panel). Nevertheless, we are still not able 
to answer the question: Which model (between these two) is the 
best to describe the data? This locates us in the context of hy- 
potheses testing problems, and we will make use of the standard 
Bayesian solution to this problem computing the Bayes factor 
for one hypotheses against the other The Bayes factor B21 for a 
model M2 against another model Mi given the data y is the ratio 
of posterior probability, namely 



B21 - 



P(y\M2) 
PiylMiY 



(1) 



the ratio of marginal likelihoods. We consider the estima- 
tion of the integrated like lihood from posterior simulations 
output JRafterv et alj l2007l) . This strategy is becoming pop- 
ular for hypotheses testing in astrophysical contexts (e.g., 
iProtassoy et al. 2002; Gregory 2005; Trotta 2008). Thus, we fol- 
lowed |^ftery_etalj (12007) and use the harmonic mean identity, 
which says that the reciprocal of the integrated likelihoods is 
equal to the posterior harmonic mean of the likelihoods. The 
simplest estimator of the harmonic mean is: 



tthmCv) 



-t — 



T-l 



(y\ff) 



(2) 



based on B draws 9^,fi,...0^ from the posterior distribution 
n{0\y). This sample might come out of a standard Markov chain 
Monte Carlo implementation, for example. Indeed, we make use 
of Monte Carlo Markov chain (MCMC) simulations to compute 
the marginal likelihoods necessary to confront one model against 
the other. 

Simulation 1: 2-Absorption-lines and Edge models without 
constraint in their parameters .— This means that, for each 
model, the MCMC routine is allowed to explore the full parame- 
ter space. We run a MCMC with B - 2000 on the Chandra data, 
using two models: 

Model 1 .- Power-law with intrinsic absorption and an ab- 
sorption Fe Edge at x 7.7 keV (quasar frame). 

Model 2 .- Power-law with intrinsic absorption and two ab- 
sorption lines at a: 8 keV and x 9.5 keV (quasar frame). 

Finally, the specification of a statistical model requires the 
form of the likelihoods terms. In our model, each likelihood term 



^/27^o-i 



exp 



(M^-yi? 



2^2 



(3) 



where v, and cri are the data (counts) and its uncertainty in the 
bin /. In this context, k - 1,2 refer evaluation of models 1 and 
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2. The index i runs over the bins falling into the range 1-2 keV 
(52 bins for the Chandra data). 

This concludes the statistical specification of our simulation. 
A major drawback of the harmonic mean estimator is its com- 
putational instability ( Rafterv et al.l 120071) . In fact, by monitor- 
ing the cumulative harmonic mean of simulation 1, we could 
see v ery large jumps, evidencing this instability (Rafte rv et alj 
|2007|) . In order to avoid statistical complications and make some 
progress, we have constrained the energies of the two lines in 
Model 2 and re-run a second simulation. 

Simulation 2: 2- Absorption-lines and Edge models with con- 
straints in the line energies .- This means that we run our 
MCMC with the energies of the two lines in Model 2 fixed at 
8 keV and 9.7 keV, values taken from best-fit parameters of C02 
and fully compatible with our own best-fit values. Here, again 
B - 2000 and the models are the same as before. 

The Monte Carlo simulations are run using parameter values 
fitted to the data under the respective models and account for un- 
certainty in these fitted values. From the resulting ;if^(l - 2 keV), 
we find that Model 2 always fits the data better than Model 1 . 
We check the stability of the harmonic means through monitor- 
ing. We note that the harmonic means for this simulation are 
stable. Finally, we compute the Bayes factor of Model 2 against 
M odel 1, u sing Eq. ^ and found B21 ~ 22 dB. |3 According 
to iJefirevsl ([1961), this can be interpreted as "strong evidence" 
for Model 2 against Model 1, given this data. In our context, 
the 2-Absorption-lines model better describes the 1-2 keV 
spectrum of Chandra, against the Edge model. 

Now, we proceed to assess the goodness-of-fit, of the two 
models presented, to the XMM-Newton data. We applied exactly 
the same methodology previously described; simulations 1 and 2 
have the same meaning, and models 1 and 2, too. But the under- 
lying data is the XMM-Newton data. The results given by sim- 
ulation 2 (this time 35 bins are included) again show Model 2 
producing smaller ;i'^(l - 2 keV), although 1) the difference is 
much smaller compared to Model 1 [only A;t'^(l - 2 keV) ~ 3]; 
and 2) the parameters are less constrained (compared with the 
same simulation using Chandra data). The computation of the 
Bayes factor this time gives B21 ~ 5 dB. In the Jeffreys' scale, 
this is interpreted as "barely worth mentioning " . This means that 
the XMM-Newton data is not able to discriminate between the 
two models. From here, it is clear that neither of the two models 
(Edge or 2-Absorption-lines model) can be applied, with 
equal success (and unambiguously) to both sets of data. In the 
next section, we will try to establish the physical scenario under 
such spectral modeling is possible, having as goal to propose 
one unified model that can reasonably present good fits on both 
spectra. 

4. Photoionization modeling of the X-ray Spectrum 
of APM 08279+5255 

We performed a photoionization modeling of the X-ray spectrum 
of APM 08279+5255, using the code XSTARQ with the atomic 
data of Bautista & Kallman (2001). The code includes all the 
relevant atomic processes (including inner shell processes) and 
computes the emissivities and optical depths of the most promi- 
nent X-ray and UV lines identified in AGN spectra. For that pur- 
pose, we built a grid of photoionization models with the col- 



^ Decibans (tenths of a power of 10), is a common unit to represent 
weights of evidence. 

'' Version 2. Ikn6. 
See |http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html| 



umn density of the ionized material {Nh), ionization parameter 
(log[^]), and Fe abundance^ as variables. Our models are based 
on spherical shells illuminated by a point-like X-ray continuum 
source. The input parameters are the source spectrum, the gas 
composition, the gas density «//, the column density and the ion- 
ization parameter. The source spectrum is described by the spec- 
tral luminosity L^ = honfi;, where Lion is the integrated lumi- 
nosity from 1 to 1000 Ryd, and J f^de - 1. The spectral 
function is taken to be a power-law f^ ~ e^", and a is the energy 
index, equal to 1 . The spectra contain edges and absorption lines 
from the following elements, H, He, C, N, O, Ne, Mg Si, S , 
Ar, Ca, and Fe. We use the abundances of lGrevesse et al.l (Il996h 
in all our models (we use the term solar for these abundances). 
We adopt a turbulent velocity of 1000 km s ', and a hydrogen 
density of 10'^ cm"^, as input parameters. 

4.1. Single-Absorber photoionization model 

Our first step was to see how well a single-absorber pho- 
toionization model can reproduce the broad-band X-ray spec- 
trum of APM 08279+5255, without any prior assumption about 
any of the parameters of interest, Nh, log(^), outflow velocity 
(vout), or abundances. So, we started by fitting a solar abundance 
model at rest in the frame of the quasar (vout - km s"'). We 
call this Model A. We apply the same photoionization model 
to both sets of data, the XMM-Newton and Chandra X-ray ob- 
servation of APM 08279+5255. The best-fit column density is 



i] X 10^^' cm-^ with log(^) 



(i.i4!;;;;^)xio^^^cm-^([i.32!0|4] 

1.50!°;°^ (log[^] = 1.50t_^-^Q% for the XMM-Newton (Chandra) 
data. The high column density is in agreement with previous 
fitted column densities to this spectrum (e.g., iHasinger et al.l 
2002.) . and with the general trend of hig h column densities ob - 
served in high-redshifted quasars (e.g.. iGallagher etal] 120021) . 
The goodness-of-fit is measured with;^'^, equal to 339.7 for 289 
d.o.f for the XMM-Newton data, and;^'^ = 241.3 (with 183 d.o.f) 
for the Chandra data. These global fits are statistically unaccept- 
able 0. Now, as we showed in the previous section, we need 
some local criterion to assess the quality of the fits in the re- 
gion 1-2 keV. Therefore, we have adopted the global param- 
eters of Model A and compute x^ (locally) from 1 to 2 keV. 



The results are 



X(i- 



2 keV) 



[XMM - Newton] = 79 for 35 d.o.f. 



and;t'^j 2 ^^^ [Chandra] =91 for 52 d.o.f. These fits are rejected 
even locally (see Table|2]for global and local ;t'^ and Pyy-). 

Motivated by the interpretation of H02 that the absorption 
feature around 7.7 keV (rest-frame) may be a Fe edge formed by 
bound-free transitions of iron ions from Fe xv to Fe xviii, and by 
the fact that the fitted edge would require Fe/O ~ 2 - 5, we have 
fitted our model to both sets of data leaving the Fe abundance 
free to vary. Again, we use a rest-frame approach (Model B). The 
result is shown in Table |2] We can see that this model provides 
a good description of the global spectrum, x^ = 310.1 (for 288 
d.o.f) for the XMM-Newton fit and;i:-2 = 204.1 (for 182 d.o.f) 
for the Chandra fit. We note a relatively important decrease of 
the column density of the ionized material, now Nh ~ 5 - 6 x 
10^^ cm"^. This is owing (in part) to the fact that the spectrum 
has strong spectral features produced by iron (in the hard X-ray 
band); for example, the edge composed by ionized species of Fe 



'' At omic abundance s are e ntered relative to solar abundances as de- 
fined in lGrevesse et al.l il996t) . with 1.0 being defined as the solar. 

^ /'289,339.7[Model A;ZMM-A'ewton] = 2 x 10"^ and 

J^i83,24i.3[Model A; Ctencfra] = 2 x 10"^, where P^.^^i is the prob- 
ability of exceeding x^ for v degrees of freedom. 
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XMM data 



Chandra data 




Table 2. XSTAR fit results and comparison between sets of data. 
XMM-Newton vs Chandra observation of APM 08279+5255. 



Observed Energy (keV) 



Observed Energy (keV) 



Fig. 2. Instrument convolved (top), and physical model (bottom) 
of Model B applied to the XMM-Newton data (left) and Chandra 
data (right). 



from Fe xv to Fe xviii, and likely Fe spectral lines. The ionization 
parameter remained almost the same. 

As we see in Table |2] Model B global-fits are fairly good 
(i.e, acceptable Py^^-^ but a closer look at the region (1 - 2) keV 
reveals important discrepancies between the model and the data. 
This difference can easily be seen in Fig. |2l In the top panels, 
we have the instrument convolved best-fit Model B (v = km 
s"') to each of the data sets, and in the bottom panels we present 
the flux spectra with the model without convolution with the in- 
strument response, so we can see physically where the strongest 
absorption features are predicted to be. In fact, there is an im- 
portant mismatch between theory and data. The statistical evi- 
dences can easily be read out from Table |2l The local j^^^j^ j^ev) 
for XMM-Newton {Chandra) is 55.5 for 35 d.o.f (65.5 for 52 
d.o.f), resulting in P(\-2) kev [Model 'Q;XMM - Newton] - 0.01 
and P{\-2) kev[Model B; Chandra] =0.1. These probabilities lo- 
cate Model B, very close to the rejection limit (a common thresh- 
old to reject the null hypotheses is Pyy. < 0.05). Based on 
the last result, we went further and explored the possibility that 
the gas absorbing X-rays in APM 08279+5255 is outflowing at 
intermediate-to-relativistic velocities. For that purpose we have 
shifted the spectra Q produced with our XSTAR-based ioniza- 
tion models, by an array of velocities, from 0.08c to 0.30c, with 
0.01c of resolution, and fit in several ways the X-ray spectrum of 
APM 08279+5255, using the XMM-Newton and Chandra data. 



Table |2] quotes the best-fit parameters of the four most in- 
teresting fits, in the context of the single-absorber model. 
Models A and B are models with Fe solar and Fe free to 
vary at Vo„, - km s"', respectively. Then, we have Model 
C and Model D, with Fe solar and Fe free (respectively) at 
Vout - 0.21c, for both sets of data, XMM-Newton and Chandra. 
The single-absorber model cannot be ruled out (instanta- 
neously), but it does not give a consistent fit to both sets of 
data (Model B for XMM-Newton and Model D for Chandra). 



^ i.e., not taking into account radiative transfer effects. 

** We made the analysis of best-fit velocity by inspecting the evalua- 
tion of ;^f2 at every outflow velocity point of our grid of velocities, while 
the other parameters of interest are varied as usual. Then, we proceed 
to adopt the velocity-model with the minimum ^2 



Parameter 


XMM data 


Chandra data 


Model A: 


^e fixed to solar; km s-' 


F 

Norm" 
Nh (cm-2) 
log(^) 

;rV(d.o.f) 

/'global [Model A] 
x2(l-2keV)/(d.o.f) 
/'(i-2keV)[Model A] 


2 J J +0.03 

1.34+jJ-«fxlO-* 
1.14t«xl023 

l-50!°:!!^ 
339.7/(289) 
2 X 10-2 
79/(35) 
4 X 10-" 


1 84+""^ 

-0 03 

1.08+o;j3x 10-" 
1.32+o"[4 X 1023 

1 50+008 

241.3/(183) 

2 X 10-3 
91/(52) 

3 X 10-" 


Model B: 


Fe free to vary; km s-' 


F 

Norm" 
Nh (cm-2) 
log(4^) 
Fe 

;r'/(d.o.f) 

Fgiobai [Model B] 
;t2(l-2keV)/(d.o.f) 
/'(1-2 keV) [Model B] 


1.32+jJ«3 X 10-" 

5.08+»-2» X 1022 
1.50+003 

5 04+^-'^^ 

310.1/(288) 

0.19 

55.5/(35) 

0.01 


J g2+0.03 

i.ii+0:«rxio-" 

5.97+038 X 1022 

1.50+°'" 

4 83+^'^^ 

204.1/(182) 

0.14 

65.5/(52) 

0.10 


Model C: Fe fixed to solar; outflow at 0.21c 


F 

Norm" 

Nh (cm-2) 

log(a 

^-/(d.o.f) 

Pgiobai [Model C] 

^2(l-2keV)/(d.o.f) 

/'(i-2kev) [Model C] 


9 19+0.03 
^•^^-0.03 

1.35!»- 2 X 10-* 

7.08!** X 1022 
1 09+0.04 

351.2/(289) 
7 X 10-3 
89/(35) 
5 X lO-** 


1.83!iJ-i!3 

i.o6+[;;j2xio-" 

6.72^4x1022 

1 79+0.08 
^ • ' "^-0.09 

246.5/(183) 

1 X 10-3 

101.2/(52) 

1 X 10-5 


Model D: Fe free to vary; outflow at 0.21c 


F 

Norm" 
Nh (cm-2) 
log(f) 
Fe 

;i'V(d.o.f) 

^global [Model D] 
;i'2(i _ 2 keV)/(d.o.f) 
/'(i-2kev) [Model D] 


2 08+""^ 

-0 03 

1.32+0 02x10-" 

3.85S X 1022 

J -76+0.06 

4 84+^-^^ 

333.7/(288) 

0.03 

71.4/(35) 

9 X 10-^ 


1 S'^+O-OJ 
l-''-'-003 

1.13+002x10-" 
3.7lto^xl022 

2 j 2+0-04 

8:36+11 

207/(182) 

0.11 
62.6/(52) 

0.17 



The error parameters are 90% confidence limits, (a) Power-law 
normalization, photons keV-'cm-2s-' at 1 keV in the observed frame. 



Therefore, we explore the possibility of a multi-component pho- 
toionization model. In Figs. |3]and|4] we plot the residuals be- 
tween each model presented here and the data. The dashed lines 
are residuals from a two-absorbers model, which we discuss 
in the next section. A detailed physical discussion of these sce- 
narios and its implications are given in Sect. |6] 



4.2. Two-Absorbers photoionization model 

Now that we have explored how the single-absorber model 
fits the data, we can go further, and see if the data sup- 
ports a multicomponent photoionization model. The simplest 
such model is a two-photoionized-absorbers model, and we in- 
vestigate if the addition of an extra component to the best- 
fit single-absorber model is statistically significant. This 
two-absorbers model consists of: one component at v = 
km S-' (rest-velocity component), and one at v = 0.18c (high- 
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XMM data 





Model D 




Model B 




Observed Energy (keV) 

Fig. 3. Residuals (in sign(data-model);tf^) of our photoioniza- 
tion models A,B,C and D to the XMM-Newton data, in the band 
(1-2) keV observed frame. Dashed line is the residual of the 
two-absorbers model. See text. 



Chandra data 




Observed Energy (keV) 

Fig. 4. Residuals (in sign(data-model);tf^) of our photoioniza- 
tion models A,B,C and D to the Chandra data, in the band 
(1-2) keV observed frame. Dashed line is the residual of the 
two-absorbers model. See text. 

velocity component) Q We take the best single-absorber 
model (from Table |2]l, selected as the best combination between 
global and local P^y.. Model B is the best for XMM-Newton 




Observed Energy (keV) 



Fig. 5. Chandra X-ray spectrum of APM 08279+5255 in the 
0.2 - 10 keV (observed frame) band. The solid thick line is the 
best-fit two-absorbers model. The inside plot shows the resid- 
ual (in sign(data-model);i'^) for the 1-2 keV band. This model 
gives good agreement on both sets of data. 




Observed Energy (keV) 



' This is the best high- velocity component we found able to fit both 
sets of data with high;i^^-probability. 



Fig. 6. XMM-Newton X-ray spectrum of APM 08279+5255 in 
the 0.2 - 10 keV (observed frame) band. The solid thick line is 
the best-fit two-absorbers model. The inside plot shows the 
residual (in sign(data-model);i'^) for the 1-2 keV band. This 
model gives good agreement on both sets of data. 



data and Model D is the best for the Chandra set. The addi- 
tion of an extra component (high-velocity component for XMM- 
Newton and rest-velocity component for Chandra) significantly 
improves the fit compared to the single-absorber model at 
the greater than 99.9% confidence level in both cases (according 
to the F-test), both globally and locally. 

First, we fit a two-absorbers model to each set of data sep- 
arately and then we do it simultaneously, to check for inconsis- 
tencies between fits. In Table [3] we have the results for these 
three schemes (Cols. 2,3,4). Figures |5] and |6] present plots of 
the best-fit two-absorbers model over the Chandra and XMM- 
Newton data, respectively. The fits between data sets give diff'er- 
ent best-fit parameters at > Icr, for six out of the seven fitted 
parameters (log ^[2] is fully consistent). The two-absorbers 
model gives acceptable fits with high ;(f^ -probability for both 
sets of data, but there exist small differences between best-fit pa- 
rameters. 

If we fit both data sets simultaneously (Col. 4), we find rea- 
sonable consistency between them and the separated fits. The 
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Table 3. Two- Absorber fits for XMM-Newton and Chandra. 



Parameter 


XMM 


Chandra 


Simult. 


XMM 


Chandra 


XMM 


Chandra 


r 


2 08+""^' 


1 84+""- 

^•°^-0.02 


2.00!-; 






2 13+"-"- 


1 01+0.02 

1.83_oo2 


Norm" 


1 41+0.03 


1 1 9+0.02 
1-1*^-0.02 


1-34!°:!!! 






1 50+0-02 

^■■'"-0.02 


1 17+0.02 

^•^'-0.02 


log^(l) 


572+0-21 

1 15+^-'^ 
3'l4+8i5 

^■^^-0.25 


5,05+0.19 

1 5r|+8:i')l 


7 34+0.18 
1 50+^-'''- 

^•-'"-0.04 










Fe 


4 70+0.23 
^■'''-0.22 


2 85+0-'o 










Nh(2)'' 


4 7i:+1.20 

t. /J Qgg 


2 74+0-68 

^■'^-1.15 


<; 4Q+0.71 

-'•^='-0.54 










log«2) 


3 03+"'5 

-'•"■-'-0.04 


3 10+0-23 


T nn+0.04 










xVid-oS) 


296.1/(286) 


186.5/(180) 


647.9/(473) 


373.8/(293) 


274.1/(187) 


298.5/(291) 


189.3/(185) 


^global [Two Abs] 


0.35 


0.38 


4 X 10-** 


7 X lO-'* 


2 X 10-5 


0.39 


0.43 


^2(l-2keV)/(d.o.f) 


44.9/(35) 


48.9/(52) 


135.1/(87) 


65.7/(35) 


69.4/(52) 


43.9/(35) 


49.3/(52) 


Ai-2k.v)[TwoAbs] 


0.14 


0.65 


4 X lO-t 


7 X 10-^ 


0.05 


0.17 


0.63 



The error parameters are 90% confidence limits. The Fe abundance is the same for both absorbers, (a) Power-law normalization, x 10 '' photons 
keV-'cm-^s-' at 1 keV in the observed frame, (b) Column density of the component, xlO^- cm-^. 



0.3 X 10 '^ ergs cm ^ s 'and 
:0.3xl0-'-^ergscm-2s-', 



most notable discrepancy is seen in the power-law component, 
with differences of ~ 10 - 15% in its parameter values. To 
check for this difference, we have taken all the parameters re- 
sulting from the simultaneously fit and fixed them to each set 
of data separately (Cols. 5,6), and compute ;^'^-probabilities. 
In both cases, the fits are rejected (in both global and local 
senses). Finally, we have taken the simultaneous best-fit pa- 
rameters and fixed them, except that we allowed the power- 
law parameters free to vary in each set of data (Cols. 7,8). 
We recover the goodness-of-fit and the model becomes accept- 
able. Then, we compute integrated observed fluxes (later in 
Sect. |6] we also compute intrinsic luminosities) using both set 
of data in the band 0.2 - 10 keV with the following results: 
F(o,2-io keY)[Chandra] = 6.9 
^(0.2-10 keY)lXMM - Newton] = 7.6 
thus the differences seen in the power-law component are re- 
flected in a change of ~ 10 % on fluxes. Pi and this could pro- 
duce small changes in the physical parameters of the absorbers 
(the most notorious are log^(l), « 30% and Nh(2), ^ 70%). 
However, the errors on the parameters allow us to obtain high 
;\f^-probability if, appart from the power-law, both observations 
are fitted with the same physical two absorbers, opening the pos- 
sibility that both absorbers have been there in both observations. 



5. Line Identification 

We are now in position of studying in more detail the possible 
identification of the feature ~ 8 keV (rest-frame) of the X-ray 
spectrum of APM 08279+5255. We will discuss two cases: i) 
the possible identification if we assume the single-absorber 
model is the best, and ii) the possible identification if we con- 
sider the two-absorbers model. 



5.1. Case(i) 

Given the good agreement between the Chandra data and our 
single-absorber model at v ~ 0.21c (Model D of Table|2l), 
we want to investigate how well constrained is the best-fit ion- 
ization parameter found and, its ability to reproduce the fea- 
ture ~ 8 keV as well as the X-ray absorption at lower ener- 
gies. For that purpose, we have taken one model in which the 



'" For completeness we have computed the observed flux of the 
~ 16 ks XMM-Newton observation (XMMl) and the result is: 



.ffO.2- 



(0.2-10 kcV) 



[XMMl] = 9.6 ± 0.6 X 10-'-' ergs cm-^ s- 



column density Nh is fixed to 5 x 10^^ cm--, built a grid in 
^, 1.5 < log(^) < 4.5 with resolution of 0.1, leaving the iron 
abundance free to vary, and computed ;t'^ produced by fitting our 
models at v = 0.21c. The main result of this experiment is that 
clearly the Chandra data strongly favors log(^) » 2.1 ±0.1. For 
instance, at log(^) a; 2.4, there exists a clear deviation of the 
best-fit model from the data with hy^ ~ 70, translated in that 
the model becomes rejected (from accepted). Thus, in Model 
D, the ionization state is mostly driven by the absorption at 
low energies, with an important contribution of the feature ~ 8 
keV. Under these physical conditions, the feature is better rep- 
resented by a complex of lines produced by L transitions of 
highly-ionized species of iron from Fe xviii to Fe xxii (Fe C2 
in Table |4]l. The centroid of this complex is located ~ 8.1 keV, 
i.e., V a 0.21(+0.01)c. 

This identification is different from the one made by C02, 
who identify this feature as Fe xxv K lines. That identification 
is not based on a photoionization model, but on searching for 
the closest (and most conservative) line (with some other phys- 
ical arguments), in likely several atomic data bases. The prob- 
lem with the identification of this feature as Fe xxv lines (in our 
Model D), is that they required that the plasma be at log(^) > 2.8. 
Only at this highly-ionized state, the ionization fraction of ions 
like Fe+^^, Fe+^^ and Fe+^^, are high enough to form spectral 
lines. In that physical scenario, a complex of lines (Fe CI), pro- 
duced by transitions of iron ions from Fe xxiii to Fe xxv, arises 
naturally, where the K lines of Fe xxv are involved. However, 
our Model D at this log(^) » 2.8 is a bad representation of the 
absorption at lower energies. In Fig.|7] we present the ionization 
fraction {q,) of Fe ions from Fe xviii to Fe xxvi of our grid of 
XSTAR photoionized clouds at different log(^). At log(^) k 2.1 
(~ best-fit ionization parameter) the predominant ions are Fe 
ions from Fe xviii to Fe xxii, and the ionization fraction of Fe+^'^ 
(?FeXxv),''5 ~ 10-^. In Table|4]we present: the transitions con- 
tributing to the complex, the laboratory wavelength, and our 
XSTAR computation of the optical depth at the line core of each 
of these lines when the plasma is at these two ionization states 
(this gives a quantitative idea of the contribution to the complex). 



5.2. Case (ii) 

Nevertheless, after proving that the addition of a second compo- 
nent to the photoionized gas, is statistically required by the fits 
in both sets of data, we are able to support the identification of 
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Table 4. XSTAR spectral lines prediction. 



Complex 



Ion 



Transition 



To(log[fl^3) To(log[g] ^ 2.8) To(log[g]«^2.1) 



FeClat 


Fexxv 


Is- 'S- Islp'F 


1.85 


110 


44 


« 10-- 


Nh = 5x\0-^ cm-2 


Fe XXIV 


ls^2s-S- ls2s2p^P 


1.86 


33 


36 


« 10-2 




Fe XXIII 


2s- 'S- ls2p'P 


1.87 


45 


108 


« 10-2 


Fe C2 at 


Fe xxn 


2s^2p- ls2s^2p^^'P 


1.88 




42 


3 


Nh = 5x10-^ cmr- 


Fe XXI 


2p^ - ls2s^2p^ ^S 


1.89 




10 


5 




Fe XX 


2s^2p^ - ls2s^2p* *P 


1.91 




3 


18 




Fe XIX 


2s^2p*-ls2s^2p^^P 


1.92 




0.4 


26 




Fe xvm 


2s^2p^- ls2s^2p''^S 


1.93 




0.1 


34 



Two different complex of lines are predicted at different ionization states. The laboratory wavelength /li^b is given in A. 




Fig. 7. XSTAR ionization balance computation (ionization frac- 
tion qf). Only shown in this plot species of Fe (xvm-xxvi). The 
single-absorber model favors log(^) ~ 2.1, and the absorp- 
tion feature at ~8 keV rest-frame would be identified as the 
complex Fe C2 in Table |4] However, the HV component of the 
two-absorbers model strongly suggest that the feature can be 
associated to the Fe xxv Kq, line (log(^) ~ 3 is precisely where 
Fe^^"^ has its peak). 



the feature as a Fe xxv line. The strongest evidence is that the 
high-velocity component is consistently fitted with an ionization 
parameter log ^(2) x 3. This is precisely where the peak of the 
ionization fraction of Fe^^'* is located. Therefore, we identify the 
absorption feature at a; 8 keV as the complex Fe CI, in which 
the major contributor is the Kq, line of Fe xxv. The physical im- 
phcations are summarized in Sect.|6] 



6. Discussion and Summary 

We have found several interesting aspects in the X-ray spectrum 
of the QSO APM 08279+5255. The physical scenarios raised af- 
te r the interpretation o f Hasinger et al. (2002) as a Fe edge, and 
of lChartas et al.l(l2002l) as relativistic (~ 0.2c-0.4c) Fe xxv lines, 
have helped to scrutinize more closely the X-ray spectrum of this 



object. If we attempt to explain the spectrum by only investigat- 
ing the full band ~ 0.2 - 10 k e V (observed frame), we find it 
hard to justify which model represents the data better in a statis- 
tical sense. After including a local criterion to evaluate the fits in 
the band 1 - 2 k e V, we were able to converge to a model that 
accounts for both low- and high- energy spectral bands. Indeed, 
by modeling the spectrum with a photoionization model, which 
gives a better (over Gaussian lines alone for instance) represen- 
tation of edges and absorption lines in the whole band ~ 0.2- 10 
keV, we are able to draw more physical information. 

We discuss the implications of our results from Sect. 
m in two parts: (i) Assuming that the X-ray spectrum of 
APM 08279+5255 can be modeled by a single-absorber 
photoionization model, (ii) Assuming that the best representa- 
tion of the spectrum is a two-absorbers model. 

(i) Changes in the physical and kinematical state of the ab- 
sorber: We see that the XMM-Newton data is represented by a 
highly-photoionized gas with column density ~ lO^"* cm"^, and 
log(^) ~ 1.5 at km s"' in the frame of the quasar. Nevertheless, 
by allowing iron to be higher than solar (model B), we obtain 
a iron abundance best-fit value of Fe = 5.04^° ^j, and the ion- 
ized gas column no longer needs to be as high as 10^^ cm"^, but 
~ 5 X 10^^ cm-^ instead, statistically improving the fit at 99% of 
significance over model A (with the same number of parameters 
than the absorbed power-law modified by an edge), similar to the 
H02 interpretation. 

But ~ 2 weeks before (in the frame of the quasar), the 
Chandra spectrum presented evidence of a material outflow- 
ing ~ 0.21c, not only showing an absorption-line-like feature at 
~ 8.1 keV (rest- frame), but also pointing to absorption at lower 
energies ~ 0.8 - 2 keV (rest-frame), coming from species of iron 
from Fe xviii to Fe xxiv, and H- and He-like ions of Ne, Mg, Si, 
S, and O. In this context, the Chandra data appears to be pointing 
to (through model D): (a) a change in the kinematical state of the 
absorbing material, or more plausible, change in the projection 
of the velocity field on our line of sight (which would yield to 
the same effect), from ~ 0.21c to ~ km s"' (since the Chandra 
observation was made first than the XMM-Newton observation); 
and (b) a slight but noticeable change in the ionization state of 
the gas, from log(^) ~ 2 to log(^) ~ 1 .5 in a timescale of ~ 2 
weeks in the frame of the quasar. However, we find it physi- 
cally hard to explain the deceleration shown by this model, in 
addition to a change in the ionization state of the gas without a 
significant change in the intrinsic luminosity of the source, and 
once we found the two-absorbers model statistically superior, 
we ruled out (i) and focused on (ii). 

(ii) Both absorbers have been there: The two-absorbers 
model consists of: one absorber at v = km s"', Nh{1) = 
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7.34+"- ^xlO^-'cm-^ log ^(1) = 1 .50+""^, and a second absorber 



at V = 0.18c, Nh(2) = 5.49+;}^^ x lO^'' cm 



-Mog^(2) = 3.00!«:«| 
(parameters coming from the simultaneous best-fit). The chem- 
ical composition of both is the same, solar in all the elements 
(see model composition in Sect. |4| except in the abundance of 
iron Fe - 2.85;^y[g. It is worth mentioning that from the sep- 
arated fits, the two observations show a change in the power- 
law of s! 10% and also show changes in the physical param- 
eters. However, they (individually) change little (i.e, no large 
difference within errors); log^(l), from ^ 1.2 to 1.5 (~ 30%); 
Nh{1) X 10^2 cm-2, from ^ 5.7 to 6.1 (~ 10%); log ^(2), from 
* 3 to 3.1 (^ 3%); Nh{2) x lO^^ ^.^^-2^ ^^^^ ^ 43 ^^ 2.7 {« 
70%, but see errors). The computation of ;^2, on the other hand, 
is very sensitive to the continuum level, and small changes in the 
power-law are easily detected by the minimization-;^^^ routine. 
We do not think this is a fundamental physical change of state 
between observations. At the moment these data do not allow us 
to discriminate whether the changes seen in the physical param- 
eters are significant (within » 2cr) or not. Apparently, apart from 
the power-law component, both observations can be represented 
by the same two-absorbers model. 

We have verified that the overabundance of Fe can be es- 
tablished respect to the abundance of oxygen (since it is the 
Fe/O ratio, which is cosmologically relevant), by comparing 
the observed spectra with the synthetic theoretical spectra at 
low-energies, if we increase the oxygen abundance. The ra- 
tio Fe/O must be a; 3 in order to obtain accept able fits in the 
low-en ergy band. O This was first noticed by 'Hasing er et alj 
(I2OO2I). and it s cosmological implications are discussed in 
iKomossa & H asinger (2003), based on chemical evolution stud- 
ies of iHamann & Ferland (,1993,) . 

This multi-component photoionization approach has been 
successfully applied to other AGNs in different b ands (X-ra 



and UV): NGC 3 783 (Netzer et al. 2003: Kronpld et al.ll2003i 
Gabel et al.ll2005l) ; Mrk 279 dFields et al.ll2007tlCostantini et alj 



2007h: NGC 4593 ISte enbrugge et al .1 l2003h: NG C 985 



(iKrongold et aP l2005h: MR 2251- 178 (TKaspi et alJ l20a 



NGC 4051 (Kraemer_etiD lIOOl lArmentrout etal] 12007 
iKrongold et al. 2007.). a survey of a sa mple of 15 Type I AGNs 
can be found in lMcKernan et al.l (l2007l) . We will treat, this multi- 
component photoionization approach, as the best solution for our 
problem. However, note that a key difference with some of these 
cited papers is that the solution for APM 08279-H5255 requires 
two unrelated absorbers at very different velocities. 

A plausible framework in which to place t hese observa- 
tional clues is in accretion disk wind models dMurrav et alj 
119951; lElvisll2000l; iProga et al. 2000; Proga 2007I) . Recent high- 
resolution spectroscopy of APM 08279-H5255, shows UV BALs 
in a wide range of velocities (5 000 < v < 12000 km s ') and 
ionization degrees (ions like C vi, O vi, N v and Si iv), im- 
plying that the outflow may be composed of multiple compo- 
nents, with velocity; densi ty and /or ionization parameter gradi- 
ents dSria nand & Petiti eanll2000l) . In particular, the model pro- 
posed by lElvisi (i2000i) predicts that the UV BAL is formed in 
the conical section of a funnel-shaped flow ranging velocities 
10000 < V < 60000 km s"', (or up to ~ 0.2c in our con- 



' ' We verified that the data is sensitive to an increase of oxygen abun- 
dance. We take the best-fit two-absorbers model and compute mod- 
els with the oxygen abundance at 1.5, 3 and 5 x solar oxygen abun- 
dance. The fit get worse with^2(0 = 1.5) = "iQA^x^iO = 3) = 333 and 
X'iO = 5) = 388. The conclusions are: 1) the data is sensitive to the 
Fe/O ratio. 2) This ratio must be « 3 + 1, in order to produce acceptable 
fits (and produce good description of the low-energy band). 



text). Furthermore, the model is very specific as to the loca- 
tion and physical properties of the X-ray absorbing material. The 
warm highly-ionized medium has column densities of ~ 10^^^ 
cmT^, and temperatures of ~ 10^ K, right pressure, and ioniza- 
tion parameter to ensure pressure equilibrium with the BELR 
clouds. Our XSTAR-based photoionization model estimates that 
the high-ionization, high-outflow velocity component (HV com- 
ponent), log(^) a: 3 has a temperature T ~ few times 10^ K, 
and the best-fit column density Nh(2) ~ lO^^'^ cm"^, nicely con- 
sistent with the main physical properties of the warm highly- 
ionized medium. The low-ionization rest-frame velocity com- 
ponent (RV component) with log(^) a 1.5 has a temperature 
T ~ 10^ K, in accordance with the t emperature of ~ few times 
10^ K found bv lKaastra etall d 19951) in NGC 5548, and the de- 
gre e of ioniz ation is consistent with log(L'^) ~ 0.3 reported by 
Mat hur et al.l (1995) for the same object, where a test to unify 
UV/X-ray absorbers in Seyfert galaxies was successfully ap- 
plied. 

Recalling the definition of the ionization parameter ^ = 
L\onl{nR^), where Lion is ionizing luminosity, n is the gas density, 
and R is the location of the absorber, we are able to compute the 
product of two un-observable quantities {nR^) from the observ- 
able ^ and Lion. We measured an intrinsic X-ray luminosity ri of 
Lx - 1.12(+0.04) X 10"*^^"' ergs s"', k being the lensing mag- 
nification factor. With this we have: {nR^)[WN]- 1.12(±0.11) x 
lO^^i cm-i and(«/?2)[RV]= 3.54(+0.11)x 10%"' cm^i.This 
is the first time this quantity is measured for this BAL QSO 
(since it is the first time ^ is properly constrained). Despite 
we still are not able to break the natural degeneracy between 
n and R (due to lack of variability), it allows us to present order- 
of-magnitude estimates on the mass-loss rate and energy bud- 
gets of the system. Assuming the HV forms part of the high- 
ionizatio n BAL flows, which re aches velocities up to ~ 0.2c (see 
Rodrigue z Hidalgo et al.ll2007l for a recent review of high ve- 
locity (HV) outflow in quasars), and that the X-ray absorber in 
APM 08279-1-5255 is pai t of th e conical section of the funnel- 
shaped flow proposed bv lElvisI (|2000) (although this estimation 
is independent of this specific geometry and can be equally ap- 
plied to a spherical shell of gas for instance, without significa- 
tively changing the final conclusions), which forms a shell with 
covering factor /cov moving at a velocity v, we estimate the mass- 
loss rate for the system: 



M„ = 4KjiH(nR)vfc, 



(4) 



where j^h is the mean mass per H atom. Using {nR^)['ii'V]- 
1.12 X lO^'*^"' cm"', an outflow velocity v = 0.18c and as- 
suming a covering factor of 0.01, we have a mass-loss rate of 
M,v a 2873^"' Mq yr ' . The total kinetic power released by the 
wind in APM 08279H-5255 is M„,vV2 = 2.6 x lO'^'^r' ergs s"'. 
The accretion rate is related with the accretion luminosity with: 



M,, 



17 



Lhoi 



0.1 



10"*^ ergs s ' 77 



■Moyr 



(5) 



where Lboi is the bolometric luminosity, 77 is the accretion ef- 
ficiency (assumed to be ~ 0.1), and c is the speed of light. 
With our measured Lx we are able to estimate Lboi - PL^- 
We adopt a bolometric correction factor fi - 20, based on 
the mean AGN SED of lEl vis etall (Jl594). Explicidy, in terms 
of the lensing magnification factor, M^cc ~ 340A:"' Mq yr"'. 



'- This is the unabsorbed intrinsic (1 - 50) keV rest-frame luminos- 
ity of the source. From the XMM-Newton data L^{XMM - Newton) = 
1.07(±0.03) X 10*' A-"' ergs s"', from the Chandra data L^(Chandra) = 
1.17(±0.05) X 10**' A:"' ergs s"'. We are taking the average of both. 
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This would lead to the conclusion that the central black hole 
in APM 08279+5255 accretes ~ 1 order of magnitude less 
matter than the estimated mass-loss rate. Recently, from high- 
resolution Chandra and XMM-Newton spectroscopic studies at 
least three systems show the property of having mass outflow 
rates ~ 1 order of magnitude larger than the accretio n rates (at 
~ 10% efficiency): the micro-quasar GRS 1915+105 ("Lee et alJ 
2002bl). fl ie Seyfert 1 galaxy MCG-6-30-15 dTurner et al. 2004: 
Sako et a l. 2003; Lee et al. 2002a), and the Seyfert 2 galaxy 



IRAS 18325-5926 (.Lee, ,2005). However, we are aware of two 
major sources of uncertainties in these estimations. Due to its 
location, the uncertainty in /cov is large because the exact ge- 
ometry in the base of the outflow is not known. On the other 
hand, the bolometric correction factor could be larger than the 
one used here, locating Mace and M^ in the same order of 
magnitude. Either we detect the outflow in a very small pe- 
riod of activity (which would be very fortuitous) or the out- 
flow plays a major role in the process of sending matter back 
to the interstellar med ium (e.g.. lCrenshaw et alj2003l:lLeel2005l: 
iKrongold et al.ll2007l) . If we assume that the mass of the black 
hole Mbh ~ 2 X 10^° Mq, (as estimated by Shields et al. 2006), 
the ratio Lboi/^edd ~ 0.8A:"' (where Ledd is the Eddington lumi- 
nosity). This number is in accordance with th e general t rend o f 
low Lboi/Leddi ratio for BALQSOs reported bv lBorosonI (|2002|) . 
With these estimations, we can derive the launching radius Haunch 
for a radiatively-driven wind: 



Haunch 



2GMb 



Ledd 



- 1 



(6) 



Here, we have used the relation vvs r for a radiatively-driven 
steady wind and assumed that at very large radii the flow has 
reached a terminal velocity v. Assuming that the wind is driven 
by a UV luminosity Luv ~ 4 x 10^^ ergs s~\ with a force mul- 
tiplier Tf ~ 100 (taken from lLaor & Brandal2002h . and using a 
flowing velocity v - 0.18c, we find Haunch ~ 7 x 10^^ cm. On the 
other hand, if we use the maximum observed velocity of the UV 
absorber (through the C vi line), for this quasar of v = 0.041c, 
we find Haunch ~ 1 X 10'** cm. One possible explanation for the 
difference in the two distances, is that the UV absorber is at a 
lower ionization state, lying in the low ionization BAL r egion 
of the funnel-shaped flowing (for visual help see Fig. 5 of lElvisI 
l200C i). and more specifically in the large radii part in which the 
UV absorber has reached its highest (terminal) velocity. A rig- 
orous study of the hydrodynamical properties of the X-rays/UV 
absorbers in AGN is beyond the scope of this paper However, 
the simplistic approach used in this analysis could be useful as a 
framework for future works on the study of the central region of 
this quasar. 

To summarize, the absorbing gas in APM 08279+5255 can 
be represented by a two-absorbers model with column densi- 
ties Nh(1) ~ 7 X 10^^ cm"2, A^/f(2) =s 6 x 10^^ cm"^, ionization 
parameters log ^(1) « 1 .5 and log ^(2) » 3, with one of them (the 
HV component) outflowing at v « 0.18(+0.01)c, carrying large 
amount of gas out of the system. The feature at ~ 8 keV (rest- 
frame) is fully predicted and reproduced by our photoionization 
model, to be a complex of Fe lines coming from high state of 
ionization, in which the main contributor is the Fe xxv /1 1.85 A, 
improving the characterization of the kinematics and the quan- 
titative evolution analysis of this high z quasar We confirm ev- 
idence for an overabundance of Fe/O, from the XMM-Newton 
observation o f this qua sar (previously inferred an d discussed in 
iHasinger et al . 2002; Komossa & Hasingerll2003L based on cal- 
cutations of Hamann & Ferland 1993). The analysis is made for 



the first time on the Chandra observation of APM 08279+5255, 
implying that both absorbers require Fe/O supersolar, placing 
similar constraints on models as before, and additionally shows 
that both independent absorbers have a similar chemical history. 
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